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Creating, manipulating and detecting spin po- 
larized carriers are the key elements of spin based 
electronics^ ^. Most practical devices^ ^ use a 
perpendicular geometry in which the spin cur- 
rents, describing the transport of spin angular 
momentum, are accompanied by charge currents. 
In recent years, new sources of pure spin cur- 
rents, i.e., without charge currents, have been 
demonstrated^ ^ and applied^^ In this paper, 
we demonstrate a conceptually new source of pure 
spin current driven by the flow of heat across a 
ferromagnetic/non- magnetic metal (FM/NM) in- 
terface. This spin current is generated because 
the Seebeck coefficient, which describes the gen- 
eration of a voltage as a result of a temperature 
gradient, is spin dependent in a ferromagnet^^^^. 
For a detailed study of this new source of spins, it 
is measured in a non-local lateral geometry. We 
developed a 3D model that describes the heat, 
charge and spin transport in this geometry which 
allows us to quantify this process^. We obtain 
a spin Seebeck coefficient for Permalloy of -3.8 
/iV/K demonstrating that thermally driven spin 
injection is a feasible alternative for electrical spin 
injection in, for example, spin transfer torque 
experiments^^. 

The interplay of spin dependent conductivity and ther- 
moelectricity was already known for half a century where 
it was used to describe the conventional Seebeck effect of 
ferromagnetic metalJ^. The discovery of the GMR effeclP 
sparked the interest of the community in spin dependent 
conductivity and no vel spin electronics which is going 
on until toda}^!^^^!!^. Due to experimental difficulties in 
controlling heat flows it was only until very recent that 
thermoelectric spintronics was investigatecP^^^ leading 
to the new fleld of spin caloritronics^^. A relevant ex- 
ample is given by Uchida et al^ who interpreted their 
results in terms of the generation of a bulk spin accumu- 
lation due to an applied temperature gradient in a ferro- 
magnet fllm. In contrast, the effect we describe in this 
paper arises from a heat current flowing through a ferro- 
magnetic/non magnetic metal junction (FM/NM) which 
creates a spin accumulation at the interface. 

The concept of how we generate a heat current over a 
FM/NM junction and subsequently measure the spin ac- 
cumulation is shown in flgure[l] The scheme is essentially 
a lateral non local spin valve structure^ with the elec- 
trical injection replaced by thermal spin injection. We 
use this non local sche me to se parate spin injection from 
possible spurious effect^^^ and because the observed 
thermally generated non-spin related voltage, which we 



refer to as the baseline resistance, allows to extract the 
temperature distribution in the device by comparing this 
to modeling^. 

We flrst formulate an appropriate diffusive transport 
theory for thermally driven spin injection. The Seebeck 
coefficient describes that an applied temperature gradi- 
ent across a conductor generates an electric fleld^^. In a 
ferromagnet, the transport processes for the majority and 
minority spin are different leading to a spin depe ndent 
conductivity a^^^ and Seebeck coefficient St,4!^^- The 
flrst is used to describe magnetoelectronics^^ in FM/NM 
systems where the latter one is usually disregarded. In 
order to consider what happens when heat is sent through 
the system, we write the spin dependent currents: 

Jt4 = -<^T4(^VMT,;-5t.;VT) (1) 

here /i^^^ is the spin dependent chemical potential. When 
a heat current Q is sent through a ferromagnet in the ab- 
sence of a charge current, a spin current Js = — = 




FIG. 1: Conceptual diagram. A charge current Jc is sent 
through ferromagnet 1 (FMi) causing Joule heating due to the 
large resistivity of FMi. The NM contacts (yellow) are highly 
thermally conductive, thereby providing heat sinks. The heat 
current Q through the center FMi/NM interface injects a spin 
current into the NM depending on the magnetization direc- 
tion Ml. The generated spins diffuse towards the FM2/NM 
interface where they generate a potential A/x = P/Hs depend- 
ing on the magnetization direction M2. As a consequence of 
Joule heating, the signal expected to arise from thermal spin 
injection scales with VT oc This potential is measured us- 
ing the indicated voltage scheme by selectively switching the 
magnetization directions Mi and M2 by a magnetic field H. 
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cri?(l — P'^)SsS/T /2 flows driven by the spin Seebeck co- 
eflicient Sq = — S^. Here P is the current polarization 
of the ferromagnet and is the conductivity of the fer- 
romagnet. 

To quantify the thermal injection of spins we consider 
the FM/NM interface and solve the Valet-Fert equatiorP 
in a fashion similar as van Son et al^ (Supplementary 
Information A). The result is depicted in flgure |2j A 
spin accumulation appears at the interface driven by the 
abrupt change of spin current going from the bulk FM 
to the bulk NM, thereby acting as an effective source of 
spins at the interface. The resulting spin accumulation 
has the following expression: 



— — ^ — = eXpSgRmis (2) 

ViFM 

where Rmis = Rn/{Rn + Rf/{^ ~ P'^)) is a conductiv- 
ity mismatch^ factor in which Ri = Xi/cii are the spin 
resistances determined by the relaxation lengths and 
the conductivities ai. For the metallic interfaces under 
consideration in this paper, this factor is very close to 
1. The resulting spin accumulation induced by the heat 
flow Q = —kpM^TFM is determined solely by the spin 
Seebeck coefficient Sg and the ferromagnetic spin relax- 
ation length Xp. Its direction is determined by the sign 



Q 



V FM 




NM 








/ A|a[ 








M 


K , 



FIG. 2: Thermal spin injection by the spin Seebeck 
coefficient across a FM/NM interface. Schematic figure 
showing the resulting spin dependent chemical potentials /i^^ 
across a FM/NM interface when a heat current Q — —kVT 
crosses it. Heat current is taken to be continuous across the 
interface leading to a discontinuity in VT. No currents are 
allowed to leave the FM, nevertheless, a spin current propor- 
tional to the spin Seebeck coefficient flows through the bulk 
FM which needs to become unpolarized in the bulk NM. This 
creates a spin imbalance /n^ — /n^ at the boundary which re- 
laxes in the FM and NM on the length scale of their respective 
spin relaxation lengths A^. A thermoelectric interface poten- 
tial A/i = P/Js also builds up. On the left side no spin current 
is allowed to leave leading to a spin accumulation of opposite 
sign. 



of the spin Seebeck coefflcient which changes sign when 
the magnetization of the ferromagnet reverses. 

Using a lock-in technique we determine the relevant 
parameters Ri(jj.V/mA) and R2{ijy / mA^) from the ob- 
served voltage^: 

V = Ril ^ R2l^ ^ ... (3) 

The baseline resistance, defined in terms of a parallel 
and antiparallel contribution as {Rf + Rf^)/2 = R\^ 
allows to extract the magnitude of Joule and Peltier 
heating effects and possible conventional Ohmic potential 
drop4^. Here R\ is determined by the Ohmic potential 
drop and Peltier heating/cooling measured by the FM2- 
NM thermocouple while the baseline resistance R2 is de- 
termined by Joule heating measured by the same thermo- 
couple. The spin dependent contribution Rf — Rf^ = Rf 
to Ri is due to a conventional spin valve signal while this 
contribution to R2 comes from thermal spin injection. 

A dedicated device was fabricated to study this effect 
and is shown in figure [3| The heating of FMi has been 
kept very localized to an area of 150 x 150 nm^ by using 
thick gold contacts. Moreover, the contacts are placed 
asymmetrically to minimize the possible current fiowing 
in and out of the FMi/NM interface. An additional con- 
tact 5 is present to be able to send a current directly 
through the FMi/NM interface. By comparing the ob- 
tained signal RI to a model (see methods), we can ex- 
tract the spin injection/detection efficiency^, which has 
been made as high as possible by keeping the size of 




FIG. 3: Coloured SEM picture of the fabricated de- 
vice. The device consists of two 15 nm thick Permalloy 
(Ni8oFe2o) ferromagnets FMi and FM2 of 1 /xm x 300 nm and 
150 X 40 nm^ separated from each other by 100 nm. They are 
connected by a 60 nm thick copper funnel with small effective 
FM/NM contact areas of 40 x 40 nm^ and 30 x 40 nm^. 5/175 
nm thick Ti/Au contacts 1 and 2 are placed asymmetrically 
on FMi to Joule heat it while contacts 3 and 4 are used to 
measure Joule heating and thermal spin injection. An addi- 
tional contact 5 is present to measure a regular non-local spin 
valve signal. 
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the FM/NM contacts small. All measurements are per- 
formed at room temperature. 

Figure [4] shows our principal results on thermal spin 
injection. Four distinct P-AP and AP-P switches are ob- 
served up to 70 nV in magnitude scaling with on a large 
background originating from the Py2/Cu thermocouple. 

The interpretation of the obtained signals requires a 
detailed knowledge of the heat, charge and spin currents 
in the device. For this purpose a 3D thermoelectric spin 
model was constructed which extends the spin-dependent 
current model^ to include thermoelectricity as well as 
thermal spin injection by the spin Seebeck coefficient. 

The calculated average contribution R2 is 2.4 /j.V/mA'^ 
slightly lower then the observed 7.69 jiiV/mA^. The dif- 
ference between the observed and modelled value was 
seen before in non-local spin valve samples^. It can be ex- 
plained by a reduction in the Permalloy thickness due to 
its oxidation, which effectively increases the Joule heat- 
ing. In the following, we scale the overall Joule heating 
in our model to fit our measured result R2. We then 
find that we were able to heat FMi to a maximum of ~ 
40K at which VT^m at the FMi/NM interface is ^ 50 
K//im. At this moment the current density is ~ 8x10^^ 
A/m^, close to the point where the device will fail due to 
electromigration. 

Electrical spin injection was also measured and the re- 
sults are shown in figure [S]^. A relatively large 9 mQ spin 
valve signal is present on top of a 1.05 background, 
being only slightly different to the 7.8 mQ and 640 mQ 
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FIG. 4: Thermal spin injection measurements, a, Mea- 
surement scheme of the experiment, b, Second harmonic mea- 
surement result i?2^^ (iiV) of the observed thermal spin signal 
as a function of I^. c, Measured second harmonic signal at 
a rms current of 1.5 mA showing the four distinct switches 
resulting from the magnetization alignment of FMi and FM2 
illustrating thermal spin injection. 



calculated signals with the metallic spin parameters Xcu 
= 350 nm, = 5 nm and Pp^y ~ 0.25 obtained from 
previously fabricated samples^ ^. Here Fpy is positive as 
shown before^^. 

The observed thermal spin injection signal = -15.6 
nV/mA^ is determined from figure |4]3. We obtain a spin 
Seebeck coefficient for Permalloy of -3.8 fiV/K^ a fraction 
of the conventional Seebeck coefficient Sp = -20 ijlV/ I^. 
This gives a polarization of the Seebeck coefficient of P^ 
= Ss/Sp = 0.19 not too different from the spin polar- 
ization of the current. It is significantly larger then the 
Sg = -2 nV/K coefficient extracted by Uchida et dlW. 
However, i t is in good agreement with recent theoretical 
prediction j-'^^l^^ l 

In addition to the thermal spin injection signal, a small 
regular spin valve signal Rf = -20 jiVt is also present and 
is shown in figure 5a. The baseline resistance Rj of 90 
liVt is in line with the calculated 95 fiVt. This is caused by 
Peltier heating and cooling of the two current injecting 
contact^. 

The negative regular spin valve signal Rf can be un- 
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FIG. 5: Rectification effects and electrical spin injec- 
tion, a, first harmonic measurement Ri for the measure- 
ment setup of figure [4] b, Calculated temperature distribu- 
tion at a height of 10 nm with a current of 2 mA sent through 
FMi. It illustrates the localized Joule heating, Peltier cooling 
and heating of the two Au/Py current injecting contacts and 
subsequent thermal conduction towards the three connected 
metallic contacts, c, Measured electrical spin injection scheme 
and resulting spin valve, d, Calculated spin accumulation at 
a height of 10 nm. A small part of the current path is short 
circuited by the Cu connection so that 4% still flows in and 
out of the Pyi/Cu contact because of its large conductivity. 
This creates a large positive and negative spin accumulation. 
Due to the asymmetry in spin injection and the asymmetrical 
placement at FM2 a small fraction of 3% is still predicted to 
give a small regular spin valve signal Rf . 



4 



derstood as follows. Due to the high conductivity of the 
copper, a fraction of the current flows into and out of 
Cu/Pyi interface electrically injecting spins. A small 
net spin accumulation at the detector interface remains 
caused by the asymmetric placement of FM2. It is il- 
lustrated by the calculation of the spin accumulation at 
the Pyi/Cu interface shown in figure |5]i which shows the 
high geometrical dependence of this effect. The observed 
is somewhat smaller than the calculated -45 /nQ. We 
believe that the small 40 x 40 nm^ size of the copper con- 
tact makes sure copper grain size, lithographic precision 
and ballistic effects start dominating the magnitude of 
this effect. 

A previous device showed a thermal spin injection sig- 
nal -5nV/mA^ at a FM-FM distance of 400 nm, only 
visible at the highest current (Supplementary informa- 
tion B). A similar calculation gives S^^ = -5 fiV/K in 
agreement with the analysis for the sample presented. 

Now that the parameters governing equation |2] are 
known we may compare this to the electrical spin in- 
jection results for the transparent Cu/Py interface. We 
can calculate that for thermal spin injection /i^/VT ^ 
2 X 10-1"^ eV m/K versus fis/J ^ ^ - 10"^^ eV mVA 
for electrical spin injection through a transparent Cu/Py 



contact. 

Due to the lateral non-local geometry and Joule heat- 
ing method used in this paper, we are limited to a maxi- 
mum temperature gradient of ~ 50 K/ jam. However, in a 
typical perpendicular geometry switching by spin trans- 
fer torque^ this does not have to be the case. In order 
to switch the magnetization by electrical spin transfer 
torque one needs a typical charge current density of ~ 
5-1011 A/m^. The same stack should be able to switch 
by applying a temperature difference of only a few t ens of 
degrees as earlier theoretical studies have indicatecP^^^. 
This simple example shows that despite the weak signals 
observed in this paper, this process can be a very viable 
alternative, or even work alongside, electrical spin injec- 
tion. 

In conclusion we have demonstrated the concept of 
thermal spin injection in a lateral metallic Py/Cu sys- 
tem. By constructing a theory for thermal spin injec- 
tion and comparing the observed potential to detailed 
3D modeling, we were able to extract a spin Seebeck co- 
efficient for Permalloy of -3.8 fiV/K. This demonstrates 
thermal spin injection is a feasible alternative for elec- 
trical spin injection in for example spin transfer torque 
experiments^^. 
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11. METHODS 

A. Fabrication 

The sample in this paper was fabricated by a 1 step 
optical and 5 step electron beam lithography process. In 
each step, metals are deposited using e-beam deposition. 
For the e-beam lithography process a PMMA 950K resist 
is used of 70-400 nm thickness depending on the thickness 
of the deposited material and resilience to Ar ion milling. 
The first e-beam lithography process produces 5/30 nm 
thick and 100 nm wide Ti/Au markers which using an 
automatic alignment procedure can be aligned to in the 
next e-beam deposition steps with high precision. In the 
next four steps, the 15 nm Py, 5/30 nm Ti/Au, 5/180 nm 
Ti/Au and 65 nm Cu layers are deposited. For the last 
three steps, Ar ion milling was used prior to deposition to 
remove any polymer residue and the Py oxide to obtain 
our highly ohmic contacts. 



B. Measurements 

The measurements were performed using a AC current 
source of a frequency < IkHz far below the characteris- 
tic thermoelectric time scale of such sized systems of ~ 
1-100 ns. The obtained signal is sent to 3 Lock-in sys- 
tems measuring the 1^^, 2^^ and 3^^^ harmonic response 
simultaneously. Care was taken in deriving Ri, R2 and 



R3 by scanning the current to make sure that even higher 
harmonics were negligible. 

C. Modeling 

We constructed a 3D model of the fabricated sample 
using the finite element program Comsol Multiphysics. 
The physics is defined in terms of a thermoelectric spin 
model where the spin up, down and heat currents are 
given by: 

Jt \ / -cr^5'^ \ / V/i^/e \ 

= - a; -cj^S^ • V/i;/e 

Q / \ -^tnt -^i^i k ) \ VT ) 

(4) 

where H^^^ are the spin dependent Peltier coefficients 
given by S^,|- Tq. Here Tq = 300K which is the reference 
temperature of the device. We take these currents to be 
continuous across boundaries. At the end of all contacts 
we set the temperature to be Tq. At contact 1 in figure [3] 
we set J^,|=J/2 to inject a charge current which is being 
sent through the system by setting /i^^|=0 at contact 2 
or 5. At all other interfaces the currents are set to 0. We 
include Valet-Fert spin relaxation by assuming VJ-^^^ = 

± '^"^~4^2^^' (^t ~ l^i) bulk. Joule heating is included 

by assuming VQ = C + ^) where a scaling factor 

C = 3.2 is used to make the model correspond to the mea- 
sured R2 . The system was meshed most accurately at the 
FM/NM interfaces where the mesh size was 1 nm in order 
to accurately calculate thermal spin injection. The de- 
pendencies R^^'^ up till R4*^ were determined by calculat- 
ing the results at ± 1 & 2 mA for the parallel and antipar- 
allel configuration. The measured resistivities cjau = 2.2 
• 10^ S/m, (jcu = 4.26 • 10^ S/m and apy = 4.32 • 10^ 
S/m were taken as inputs for the model. In this model, 
the substrate was also taken into account^. The Seebeck 
coefficients 5'a^=1.7 /iV/K, 6'cn=1.6 /iV/K, Spy = -20 
/iV/K and thermal conductances k^^^ = 300 W/m/K, 
kcu = 300 W/m/K, kpy = 30 W/m/K, ksubstrate = 1 
W/m/K were taken from various sources in literatur d^ l ^^ l 

III. SUPPLEMENTARY INFORMATION A 

Here we calculate in more detail what happens when 
heat is sent through the FM/NM system in figure |6] We 
begin by writing the spin dependent currents: 

J^^l -cr^,;(^V/it,; - 5'-^,;VT) (5) 

here /i^^^ is the spin dependent chemical potential. When 
a heat current Q is send through a bulk ferromagnet in 
the absence of a charge current, a spin current Js = 
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— J| = apil — P^)S'sVT/2 flows, driven by the spin 
Seebeck coefficient Ss = — S^. Here P is the current 
polarization of the FM and ctf is the conductivity of the 
ferromagnet. Charge and spin current conservation^ ^ 
leads to the thermoelectric spin diffusion equation: 

VV. = g+e(^(VT)2 + 5,V2T) (6) 

where /ig is the spin accumulation /i^ — /i^. In addition 
to the Valet-Fert spin diffusion equation V^/ig = two 
source terms are present. Both terms can create (albeit 
small) bulk spin accumulations. 

In figure [6] we send a heat current Q through the 
FM/NM interface while we allow no charge or spin cur- 
rent to leave. The heat current Q = —kVT needs to be 
continuous throughout the system, leading to VTpM = 
kNM /kpM^TNM at the interface. Since VT is constant 
in both regions individually, and for first order effects we 
may assume Ss is constant, the source terms in equation 
[6] are irrelevant. Therefore, we may use the standard 
Valet-Fert spin diffusion equation to solve the bulk spin 
accumulation leading to the general expression for the 
spin dependent potentials in the bulk: 



with A-D the parameters to be solved in both regions. 
At the FM/NM interface we take the chemical poten- 
tials /i^^i to be continuous as well as the spin dependent 
currents J^^^. At the outer interfaces we set the spin de- 
pendent currents to zero. This leads to a set of equation 
which can be solved. We obtain: 

B = e^i^^i^VT^M ^ cSfmS/Tfm (8) 

where we use the definition of the conventional Seebeck 
coefficient of a ferromagnet Sfm^- The spin accumula- 
tion at the interface is: 

— = eXpSgRmis (9) 
y iFM 

where Rmis = Rn / {Rn + Rf /(I ~ ^^)) is a conductivity 
mismatch^ factor in which Ri = Xi/ di are the spin resis- 
tances determined by the relaxation lengths and the 
conductivities cr^. 

IV. SUPPLEMENTARY INFORMATION B 
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FIG. 6: Thermal spin injection by the spin Seebeck 
coefficient across a FM/NM interface. Schematic figure 
showing the resulting spin dependent chemical potentials /i^,^. 
across a FM/NM interface when a heat current Q — —kVT 
crosses it. Heat current is taken to be continuous across the 
interface leading to a discontinuity in VT. No currents are 
allowed to leave the FM, nevertheless, a spin current propor- 
tional to the spin Seebeck coefficient flows through the bulk 
FM which needs to become unpolarized in the bulk NM. This 
injects a spin imbalance /i^ — /i^. at the boundary which relaxes 
in the FM and NM with the length scale of their respective 
spin relaxation lengths A^. A ther moelectric interface poten- 
tial A/i = P/Js also builds up(^. On the left side no spin 
current is allowed to leave leading to an opposite spin accu- 
mulation. 



A SEM picture of the previous sample is shown in fig- 
ure [t] (a). A regular spin valve signal was measured by 
sending a current from contact 1 to 3 and measuring the 
potential between contact 5 and 4 of which the result is 
shown in figure l7|(b). In this case a 13.8 mQ background 




B (mT) B (mT) 



FIG. 7: Previous device results, a) Coloured SEM figure 
of the device. The sample consist of the same two ferromag- 
nets which are now placed 400 nm apart. It is connected by 
a copper V shape instead of a funnel, b) Non local spin valve 
signal by sending current from contact 1 to 3 and measuring 
the potential from contact 5 to 4. c,d) Thermal spin injection 
result. The current is now send from contact 1 to 2 while the 
potential was measured between contacts 5 and 4. 
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Ri is observed on top of a non local spin valve signal 
Rf of 3 mQ. The background is originating from Peltier 
heating/cooling of the FM/NM interfaces^. Both signals 
are close to the calculated 14.1 mQ and 4.1 mQ. 

When the current is send from contact 1 to 2, we ob- 
tain the results shown in figure [t] (c,d). A regular spin 
valve signal Rf of 10 jaQ is observed on top of a small 
-15 background R^. This is somewhat different then 
the calculated -100 fiQ. background and -4 fiQ. However, 
these effects are highly dependent on the exact geometry 
and are due to the small 30 x 30 nm^ size of the con- 
tact. This makes sure grain size, lithographic precision 
and ballistic effects dominate. 

Thermal spin injection was observed and is shown in 
figure [7](d). The background R^ is again larger then the 
calculated 3.4 /j.V/mA'^. If we compensate for this in the 
modelling we obtain from the observed ^ -7 nV/mA^ 



signal a spin Seebeck coefficient for Permalloy of ~ -5 
I^V/K. 

We note that our thermal spin injection signals i^^^^ 
cannot be due to an induced temperature dependence 
of Rll because R^I^, representing the Joule heating in- 
duced change of Rfl should be ^10 times larger then 
the possible Peltier heating induced i?!^^ at the currents 
we used in both samples. However, Rsl^ was simulta- 
neously measured and found absent. Furthermore, The 
signal Rfl of both devices is of different sign, excluding 
any further dependence of i^l^^ on Rfl. 

We conclude that also in this device we have good 
agreement between observed and calculated thermoelec- 
tric voltages when we apply a similar correction for the 
Joule heating. A very similar value for the spin Seebeck 
coefficient was found. 
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